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O ' Abstract 

CN ! The Bogoliubov procedure in quantum field theory is used to describe a relativistic 

o 



almost ideal Bose gas at zero temperature. Special attention is given to the study 
^ ■ of a vortex. The radius of the vortex in the field description is compared to that 

I . obtained in the relativistic fluid approximation. The Kelvin waves are studied and, 

^Jjl for long wavelengths, the dispersion relation is obtained by an asymptotic matching 

method and compared with the non relativistic result. 
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1 Introduction 

An interesting issue in the study of relativistic vortices involves the similarities and differ- 
ences between global strings and relativistic superfluid vortices [H I2l 01 IH IS] • 

This paper deals with a relativistic treatment of the almost ideal Bose gas at zero 
temperature and its approximation as a relativistic fluid. More specifically we study a 
cylindrical vortex and its Kelvin waves. 

We know that an almost ideal Bose gas at zero temperature is well described by the 
Gross-Pitaevskii equation [HI 13 El which governs the evolution of the macroscopic wave 
function caracterising the condensate. This suggests (Sec. 2) to start with a quantum field 
whose potential m^0*0+A(0*0)^ does not exhibit spontaneous symmetry breaking, in order 
to represent weakly interacting relativistic Bose particles. By the Bogoliubov procedure we 
show that the classical field equation can be interpreted as a "relativistic Gross-Pitaevskii 
equation". In order to apply the Bogoliubov procedure we have to verify that the mode. 
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in which the bosons condense, reahzes the minimum of energy with the constraint that the 
conserved charge is fixed by the number of particles. This is done in Sec. 3. 

This prehminary study makes things very similar to the non relativistic counterpart 
(cf. for example and allows us to interpret the classical field equation as representing 
the motion of the condensate of an almost ideal gas at zero temperature. 

In Sec. 4 we establish the approximation of the field equations by a perfect isentropic 
irrotational fiuid in the formulation of Lichnerowicz JU]. In Sec. 5 the cylindrical stationary 
solution(that is a vortex) is studied. We begin with the exact solution in the approximate 
fiuid theory and compare with an asymptotic solution of the exact theory. 

In Sec. 6 we study the Kelvin waves. The oscillations of a rectilinear vortex were initially 
studied by Lord Kelvin [11^ in the context of the classical fiuid. In Bose condensate theory 
the oscillations of a non relativistic quantum vortex were first discussed by Pitaevskii [7j 
who found in the long wavelength limit the dispersion relation: 

fl2 1 

a; = — In (1) 

2m qro 

where tq is the core radius. This result is refined by Grant who has found: 

2 

uj = —{ln 7-0.115} (2) 

2m qro 

where 7 is the Euler constant. Davis and Shellard jl] have shown that a spinning global 
string behaves like a vortex in a superfiuid medium. The dispersion relation of the Kelvin 
modes of such a vortex was given by Ben-Ya'acov In the non relativistic limit, he 

obtains 

u; = |^{lnl-7 + 1.5} (3) 
2m qe 

where e is a cutoff parameter which depends on the core radius parameter S of the string 
by 

In 4 = 1.615 -In 2. (4) 



This last formula is obtained for a vortex ring. Comparison of (j21) and Q yields exactly 
the same formula if tq = 6. 

In this paper we cannot follow the method of ^2] which uses explicitly an equation of 
motion of the vortices based on the spontaneous breaking symmetry of the potential. We 
shall rather transpose the method of the matched-asymptotic-expansions of Roberts and 
Grant [HI [T2] to the relativistic field equations. The result is 

q2 



^-Vl--o^|;^il---7-0.115} (5) 

where vo is the velocity at r = tq of the classical relativistic fiuid associated with the 
condensed bosons. Finally in Sec. 7 we discuss some of our results, in particular we compare 
((S)) with the dispersion relation of the Kelvin modes of a global string. 
In this article the signature of the metric is ( — h ++) and c = h = 1. 
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2 Classical field equation and Bogoliubov procedure 



Let us consider a phenomenological system of electrically neutral bosons described by a 
complex scalar field. It is worth noting that a complex field (non Hermitian operator) can 
describe particles and antiparticles without electric charge (cf. Landau and Lifchitz [T3]). 
In this case the conserved charge Q is associated with the difference between the numbers 
of particles and antiparticles. 
The lagrangian density is: 

C = -9^0*5^0 - V{(P*(P) (6) 

with potential 



l/(0» =mV</> + A(0*0)^ (7) 

The constant m? is positive and represents the squared mass of the particles, the coupling 
constant A is also positive so that the interaction is repulsive. If A is small we can describe 
an almost ideal relativistic gas of bosons at temperature T = by the field 

^ = + V^' (8) 

where is a classical field which will represent the condensate; it verifies the Euler-Lagrange 
equation 

- V 



- d^d^cf) + m20 + 2A0*00 = (9) 



and tp' is a fluctuation. 

If we suppose that the system is enclosed in a box of volume V , the field can be expanded 
in a Fourier series 



px = —p^t + p • X 



with 

and 

e = p^ = ^p2^^.^ (11) 

where eo is a constant whose precise value will be fixed below. The operators Cp, bp and 
their conjugates CLp,bp are the annihilation and creation operators for the particles a and 
antiparticles b: 



;] = h KK] = ^- (12) 

Let us note that the Fourier expansion (jlOj) and the commutation relations (fT^ usually 
given for free fields are valid also for interacting fields (cf. Lee [T^) . Of course in this 
case Op and bp and conjugates can have a complicated time dependence in contrast to free 
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fields where these operators are time independant if we choose for eo the mass m of the 
free particle. 

For a uniform weakly non ideal gas in a large volume V at temperature T = 0, almost 
all the N particles, A^^o ~ N, condense in the zero mode of the particles. As it will be 
justified in the following section, this mode represents a minimum of energy when we 
impose the constraint that the conserved charge Q is fixed: Q = Nq. Hence, since A^^o 
is a macroscopic number, following the Bogoliubov approximation ^7|, we can negect the 
commutator [oq, a^] = 1 as compared with the huge eigenvalue A^'o of ag Oq, so the operators 
ao, ag can be approximated by classical numbers 

ao ^ a+ ^ V^iVo. (13) 
Therefore the Fourier expansion (jlOj) can be decomposed 

^ = ^=^e-o* + ^' (14) 



where the classical field 



^" e-»* (15) 



represents the condensate of the particles of zero momentum and the operator ip' referring 
to the remaining non condensed particles is considered as a small perturbation. 

Let us note that we have not fixed the arbitrary constant eo. Since 0o given by ()15|) 
must be a classical solution of the equation (jHI), eo, supposed positive, is fixed by 

eo(eg-m^) 

no = (16) 

where 

^o = Y (17) 

is the constant density. 

Let us remark that we consider a gas of bosons at low energy. More precisely it is a 
relativistic gas but a weakly relativistic gas where creation and destruction of antibosons 
are negligible. Of course we could also consider the symmetrical situation of a weakly 
relativistic gas of antibosons living in an antimatter world. In this case the constraint is 
Q = —Nq where A^o is the number of antibosons, the condensation would be on the zero 
mode of the antiparticles and the classical field representing the antiparticle condensate 
would be given by (fT3j) with opposite phase. 

The generalisation of the Bogoliubov approximation to the case of non uniform and 
time dependent configurations is given by (jH)) where the classical field 0, which represents 
the condensate, verifies the equation Q. This completes the comparison with the non 
relativistic theory (cf. for example jUj). 
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3 Minimum of energy at fixed conserved charge 

In this section we show that the homogeneous solution (fTH|l and (fTTHl reahses the minimum 
of energy when the number of particles is fixed. 

It will be convenient to write the classical complex field in polar form: 

0(x) = F(x)e^^(^). (18) 

The Lagrangian (jH]) becomes 

C = -d^Fd^'F - F^d^Sd^'S - V{F^) (19) 
and the Euler-Lagrange equations take the form 

- Z^iZ!:^ + d^Sd'^S + V\F^) = (20) 
F 

where V'{F'^) is the derivative 

V^iF^d^S) = 0. (21) 
The 4-current and the energy density are given respectively by 

J, = -i{<j)*d^<j) - d^<P*<P) = 2F'd,S (22) 

and 

n = {doFy + F^idoSy + (VF)2 + F^VSy + V{F^). (23) 
Let us introduce the conjugate momenta of the fields 

dC 

nF = —r = 2doF (24) 

and 

BC 

^s = j^= 2F%S. (25) 
The charge and the Hamiltonian can be respectively written 

Q= f fd^x = [ -2F^doSd^x = [ -nsd^x (26) 
Jv Jv Jv 

and 

H= (-4 + 2^ + (VF)2 + F\VSY + V{F^))d''x. (27) 

The extremum of the energy H, when the conserved charge Q is fixed by the constant 
Nq, is determined by the variation 

6{H-i2{Q-No)) = (28) 



where /x is a Lagrange multiplier. The variation on irp, its, F, S, gives respectively 



= 0, (29) 

1 TTq 

2;^ + /^ = o, (30) 

that is 

F = F(x), (31) 

S = -fit + So{^), (32) 

and implies that these functions are solutions of the equations (f^ . (jSIJ. 

Bekenstein and Guendelman show that this type of solution is a local energy 
minimum within its sector of fixed charge. In particular the homogeneous solution (fT^ 
and (fT^ belongs to this class of solutions with fi = eo and realises the minimum of energy 
when the conserved charge is fixed by the number A^q of particles, that is 

3 1 

H = ^iVoeo + -No—. (33) 

4 4 eo 

We can show that (jH^ is effectively the energy minimum as follows. A solution of 
can be written as a perturbation of the homogeneous solution (fT3j) . (fTBj) : 



^ + /(x), (34) 

S = -eot + 5o(x) (35) 

where /(x) is constrained by 

J^(2^J{^) + f{^)y'x = (36) 

in order to maintain Q = Nq. 

If we introduce f|M| and in several terms are canceled by using and we 
obtain finally 

rr 3,, 1.. 
= ^iVoeo + -^No- 

+ I (^(V/(x))2 + (y^+ /(x))^(VSo(x))2 + A (2y^/(x) + rf'^37) 

The integral on the right hand side is obviously positive therefore ()33p is the minimum 
with Q = No. 
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4 Fluid description 

We shall express the modulus F{x) of the field in the form 

Below it will turn out that n{x) is the proper density of the particles and h{x) the enthalpy 
per particle in the fiuid approximation. The expression (jHHjl can be suggested by the 
homogeneous solution (jlSp where ^ = is the constant density and eo an energy scale 
fixed by the relation |T6|l . 

The 4-gradient of the phase d^S can be written as a covector in terms of a unit 
4- vector u^: 

= d^S = hu^ (39) 
hence from ()38|) the 4-current ()22|) takes the hydrodynamics form: 



= nu^ (40) 

where we interpret n (if it is positive) as the proper density and as the velocity 4- vector 
of the fiuid. 

From we deduce also 

d^C, - d,C^ = (41) 

and 

C^C^ = -h^ (42) 
Equation ()21|) is the familiar continuity equation 

y^.3^ = 0. (43) 



______ ^ ^ I-'' F 

In equation ()2Up the term — can be neglected if F varies slowly compared to 5*. This 
term is analogous to the "quantum pressure" ^ ^^^ in the non relativistic theory [0]. 
Let us look at this approximation more closely. The equation fl20|) takes the form 

df.Sd^'S + V'{F^) = 0, (44) 

hence we have 

71 

- Ct'C^ = m^ + \- (45) 

which combined with (j42|) gives 

A 

This relation n{h) is identical with the relation no(eo) given by but (j46j) is valid within 
the approximation considered here. 
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The equations (PT|l . (P|l . (P|l and 



= -^^M (47) 

are the relativistic equations of a perfect isentropic irrotational fluid ( cf. JUl and 
where these equations appear in the same notations). For the most general framework in 
relativistic hydrodynamics (multicomponent fluids, superfluids) , the reader is refered to 
[20] and [21].) 

The nature of h is now clear. It is the enthalpy per particle (with the mass energy m 
included). The equation is the equation of state which is automatically imposed. 

If the quantities p and p designate respectively the proper energy density and the 
pressure, the enthalpy per particle can be written 

h = P±^. (48) 
n 

Since the entropy is zero we have 

dp = ndh (49) 

The equations pH|l and define p and p in function of h and n. combined with 
can be integrated: 

P= ^^'-^f^' =\F\x). (50) 
The velocity of sound Vs = \/ dpjdp is given by 

n(h) —w? 1 , 



We remark that — > 1/3 when m ^ which is the the correct ultrarelativistic limit. 

The expressions fl48|) and (jSUjl can also be obtained by identifying the energy-momentum 
tensor of ()19p (in which the derivatives dfj_F are neglected), that is 

T^, = 2F%Sd,S + r^,,{-F^dpSdPS - V{F^)) (52) 

with the perfect fluid energy-momentum tensor. Finally we notice that the approximate 
equations can be obtained from a generic Lagrangian proposed by Carter to describe 
the irrotational fluid in a condensate: 

£ = -^F^d^Sd^S - U{F). (53) 

In our case: 

U{F) = \{m''F^ + \F^). (54) 
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5 Cylindrical vortex 



5.1 Vortex solution of the fluid approximation 

It is convenient to begin with the approximate theory governed by the equations (j^T| . (02j, 
(gSl) and (jUj) with the equation of state of the fluid (gHl). 

We know (cf. [TH], [23], that in cyhndrical coordinates 

ds^ = -df + dz^ + dr^ + r^de^ (55) 
a stationary cyhndrical solution is given by 

C^ = (-E,L,0,M) (56) 

where E, L, M are the three constants associated with the three symmetries. By a Lorentz 
transformation along z we can always choose L = and have a circular flow. 
The enthalpy h is a function of r given by 

= -C^C, = E^-L^-^. (57) 

The equation of state of the fluid (|46|) yields two values h = and h = m for which 
n vanishes. When we move from the exterior toward the center of the vortex, h{r) is 
decreasing; h is always larger than m, so n vanishes first at h = m which determines the 
radius 

M 



m 



2 



of the core of the vortex. Below this value n becomes negative, which is meaningless in 
this approximation of classical fluid. Let us observe that with a familiar polytropic fluid 
the core is empty, n is zero. The sound velocity vanishes also on the radius vq. 

We shall now suppose for simplicity that we have a circular flow with L = 0. Using 
PUj) . (jSni), (j37j) . we write the wave function of the vortex of the approximate equations: 



2A 



,i{-Et+Me) 



(5J 



where we have suppressed an undetermined constant phase. 

To make this wave function single-valued, M must be an integer, noted v. The domain 
of validity of (jHH|l is r > ro where 

ro = ,J , (59) 



m 



is the core radius for u = 1. The length ro sets the limit below which we cannot neglect 
the effect of the quantum term. 
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It is immediate to write the velocity of the circular flow 



V 



Let Vq be this velocity at r = Tq for t/ = 1, we have: 



E = -j^=. (60) 



5.2 Vortex solution of the exact theory 

A vortex solution of the exact equations and (PT|) must verify (jUJ, (jl^ . (jlHj) and 
(I47|) whose stationary cylindrical solution is always given by (|56j) and (|57p. The relation 
between n and h is no longer (pUj) but is a function of r and = is also a function 
of r. So the vortex = Fexp iS is determined (with L = 0) by 

S = -Et + uO (61) 

and by the the differential equation for F{r) deduced from ()2()|1 : 

1 

-dr{rdrF{r)) - -^F{r) + {E^ - m^)F{r) - 2AF3(r) = (62) 

This type of equation has been discussed in the literature many times (see for example 
|H]). Analytical solutions are not kown, but it is easy to discuss its asymptotic properties 
for r — >■ and r — > cxd and to compare with the solution of the fluid approximation. The 
presence of the term z/^F/r^ means that F{r) must tend to zero when r — if the energy 
in the volume V has to be finite. For small r the cubic term 2XF^ is the smallest and can 
be neglected, the solution is the Bessel function 

F{r)=AJ4{E^-m^y^^r] (63) 

as r — i> 0. 

As r oo the first two terms of (j62j) cancel: 

F{r) =Fo = (64) 

An approximation to the solution is obtained by matching the functions ()63p and ()64|) 
and their derivatives. Let r = oq be the radius of the junction. The constant A is fixed by 
the continuity of F at r = (we are not concerned by this value) and by the continuity 
of the derivative ^(oq): 

Jl[{E^-m^f/'a,] = Q. (65) 
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So Oo is given in term of the first zero of J^, let say Zy. For the winding number z/ = 1, 
Z\ = 1.84. The radius Oq of the junction, which is naturally interpreted as the radius of 
the vortex, is given by 

1.84 , , 

We observe that tq the radius of the core of the vortex in the fluid approximation and 
Oq differ only in a numerical factor close to unity, they are typically of the same order of 
magnitude. The length scale oq ~ tq is of the order of the so called "healing length" jUj. 
This result indicates that the fluid approximation is a good practical model for the study 
of the vortices in a relativistic Bose gas condensate. 

6 Kelvin waves on the cylindrical vortex 

We consider now a generic perturbation of the vortex (p = F exp {iS) determined by (|?)T|l 
and dni: 

= + A0~e^^(F + AF + iFA5). (67) 

Substituting ()(i7|l in the Euler-Lagrange equations and assuming that the perturbation is 
small, we can linearise in AF andin FAS. Then it is possible to expand AF and FAS' in 
wave modes and to consider each mode separately: 

AF = M{r)cos{-ujt + qz + 9e), (68) 

FAS = N{r)sm{-tut + qz + ue). (69) 
M(r) and N{r) verify the linear equations 

M" + -M' + {E^ + u^-q^- 6XF^ - m'^)M - ^[{u'^ + v'^)M + 2uvN] + 2EujN = 0, (70) 

N" + -N' + {E^ + uj^-q^- 2XF^ - m^)N - \Uiy^ + u^)N + 2uvM] + 2EujM = 0, (71) 
We can also |13^ introduce equivalently 
AF + iFAS = H+{r) expi{-ujt + qz + v9) + H^{r) exp -i{-ujt + qz + vO) (72) 

with 

which verify the linear equations: 

-{rH'J + [{E ± ujf -q^- = (m^ + A\F^)H^ + 2\F^H^ (74) 
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Since we are interested in the small motion of the central line of the vortex, let us concetrate 
in the neighbourhood of the origine. For small r, F{r) is given by (jUHj) hence 

F{r) ~ Ar\ (75) 

We can neglect the two last termes of the r.h.s. of Eqs. ()74|1 which are approximed by Bessel 
equations whose the solutions are 

= 5J|,±p|[((E±cu)'-g2-m2)5r] ocrl'^^^l (76) 

The center of the perturbed vortex is determined by the zero of the solution and since 
F{r) vanishes at the origin, Eq.()67p requires that at least one of the functions H± is 
different of zero at the origin. Suppose it is we must have 

u + u = Q (77) 

and consequently if_ oc r'^'^ is negligible. Hence near the origin we get 

= expi{-Et + vO) {Ar" + H+{r) expi{-ujt + qz + vO)) . (78) 

Let the constant B = —ay A, eq. (fTSjl becomes 

= A exp i{-Et) [r^e'"'^ - a" expi{-ujt + qz)) . (79) 

In the case u = 1 the cancellation of is obtained for 

r = a , e = -ut + qz (80) 

which is the helical motion of the center of the vortex, that is the Kelvin waves. 

,2 



Let us return to the equations ()70p and ()71|). If = and if we change m? into — m' 



they are identical to the corresponding equations for the global strings of [3^ . We can verify 
that for z/=l,z/ = — l,a; = and g = there is an exact solution: 

Fir) 

M{r) = F'{r) , N{r) ■- 



r 

In the following we shall have always z/ = — z/ = 1 since we are interested in Kelvin waves 
and we shall use non dimensional variables. 
Let 

y = -, (82) 

F{r) = Fo/(-) = Fofiy) (83) 

where ro and Fq are given respectively by and (jMj) . 
The equation (jHH) becomes the standard equation 



r + -/' + (1 - -f = 0. (84) 



1 1_ 
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For large y the solution is 



1 9 161 
/(2/) = l-7U^-^-T7^ + ---' (85) 



2j/2 8i/4 IQy 

for small y we have 



f{y) = ky + 0{y') (86) 

where is a constant which can be obtained numerically. 
From the equations and (jUUI) we have naturally 

M(r) = FoM(i/) and N{r) = FoN{y). (87) 

We introduce also lj = uvq, q = qvQ, E = Evq, rh = mro, t = rot, z = tqz. The equations 
(f7n|) and (dH) becomes: 

M" + -M' + (cu^ - g2 + 1 - 3f{y))M - ^(M - N) + 2EujN = 0, (88) 

y y 



N" + -N' + (uj'^ + f(y))N - —(N - M) + 2EljM = (89) 

y y 

Since u = — z/ = 1 the solution of these equations can express the deplacement of the 
vortex. This is the case of the exact solution: 

Mo = ny) , No = ^ (90) 

y 

obtained for u = q = 0. 

We can now examine the solutions of Eqs. (jHHj) and (|H^ and compute the frequency lj 
in the limit g — of the long wavelength oscillations. To solve the equations we shall follow 
the method of "matched asymptotic expansions" used by Roberts and Grant [HI [T2] . 

We consider first the solution in the inner cylindrical region of radius y centred on the 
oz axis in the form of an expansion about the exact solution (jMIj) . This solution is then 
examined for y ^ oo. We consider also the solution in the outer region caracterized by 
a new stretched coordinate s = qy which is of order unity for large y. This solution is 
examined in the limit s 0. Finally the two asymptotic solutions (y — » oo and s — 0) 
are matched in an overlap domain, where both the inner and outer expansion are valid. 

From this matching we shall obtain the dispersion relation of the Kelvin waves. 

6.1 Interior solution 

We expand M and N around the exact solution Mq and Nq: 

M = Mo + fMi + ■■■ (91) 
N = No + fNi + ■■■. (92) 
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Then, substituting and into the equations (|HHj) and (|H^ we obtain, by equating 
the coefficients of 



+ -M[+ (l-^- 3fiy)) Ml + ^iVi = f'{y) - 2Euj, 

y \ y J y 



m 

y 



(93) 



N'l + -N[+(l--^- fiy)] iVi + = ^ - 2Eu;^f\y), (94) 

y \ y^ J y^ y 



with 



u = oJiq^. (95) 

To solve and (jMj) we substitute /(y) by the expansion vahd for large y. 

First we search for solutions for the homogeneous system which is transformed in a 
linear equation of 4*^* order in Mi. 

.-Ai) 6 -(3) , „ 7 8 76 , 10 1 16 380 

+ -Mf + -2 + — + — + — + ■■ ■ M' + + ^ - ^ ^ + ■ ■ ■ M 

y y^ y^ y^' ' ^ ' y y^ y^ y^ 

6 9 18 

+ (-^-4-4 + ---)Mi = (96) 

y y y^ 

Infinity is an irregular singular point of rank of this equation. The roots of the charac- 
teristic equation are ±-\/2 and 0. The root has a multiplicity 2. By applying standard 
techniques about asymptotic solutions of linear equations [25j, it is possible to obtain 
asymptotic expansion of four independant solutions: 



n — 1 I 9 I 483 I 
~ 1^3 2y5 Sy"^ 

_ 1 I 21ny I im_ I 9]ny , 2]_ , 483 In^; , 
^"^ ~ y 36y3 ?/5 "^4^5"^ 4j^7 ■ ■ ■ 

C/3 - exp y/Zy)^yi -h 255^2 6144y2?;3 + 3932162/4 + 

^ _ ovT^ 1 n -U _5 L 65 I 5981 I 181825 , 

From these solutions Mi we deduce the corresponding solutions Ni. 



(97) 



f ^1 1 9_, 

•'I y 2y^ 8y5 ~'~ 

f = 1,-L ihl _L ^ _ iM ZZ ilM _L 

J 2 y "T -r op 3 10 3 . 5 -r 



y 2|/3 8y5 
2 In; 

y ' 36y 18y'^ iy^ 

h = exp (-v^y) A.(_ 1 + -f3 - + ■ ■ 



F ^ 8V2;/3 25%5 

/4 = exp(v/2,)^(-^-^-^ + 



(9J 



From the solutions (j97|) and (j98j) we can construct the fundamental matrix of the 
homogeneous system. Then by the variation-of-constant formula, we obtain the solution 
of the inhomogeneous system and This solution depends on four arbitrary 

constants. By choosing two of them we can suppress in the solution the exponentials 
exp (±^/2y) which are not compatible with the exterior solution. Finally, we obtain: 

X In 1 

Ml = -— + -{A + Eul) + ■■■, (99) 

2 y y 
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iVi = Iny + A,y + 1 In^ y + 1^ (2A, + 2Eu, + ^\+^ + ... (100) 
2 2?/ ?/ V 4:J y 

where Ai and are the two remaining constants which can be determined by matching 
the interior solution and the exterior solution. 

By combining the equations and ()94|1 we can obtain an integral relation which will 
be useful later in determining u. If we multiply by yf'iy) and by /(y) and add 
the two resulting expressions, using Eq. (jH^ we obtain: 

^ i-fNi + fN[ - f'N, + yf'M[ - yfM, + yfM, + f'M, - -fM,] = 

dy \y y J 

-2Eu,{f)' + y{ff + -f (101) 

y 

We integrate from to 00 and we obtain after substituting the value of Mi, Ni and / for 
y ^ 00 and y — 0: 



2Ai + ^ = -2EUJ, + y iff dy + jim (^^^ ^fdy - In A 



A 

(102) 



In obtaining this relation. In 00 has been canceled from the r.h.s. and the l.h.s. If ^ = 1 
the relation p02j) is identical to the corresponding non relativistic relation of Grant 
The numerical value of the integral part is known: 



y iff dy + Jim (^^^ ^fdy - In | = -0.115 



(103) 



6.2 Exterior solution 



Let us return to the equations (jHH|) and (jH^ and expand them for large y. Replacing y by 



I where s is a streched coordinate since q is small, the equations become: 



_2d^M q^dM 
^ ds^ s ds 



^ir - r + 1 - 3 1 



S2 



M 



2^ {M-N)+ 2EuifN = 0, (104) 



2 



S2 



+ 



2-4 -2,1 / 1 

^ig -g +1- ( 1-^ 



^ d'N ^ q'dN 
ds'^ s ds 

-2\ {N-M)+ 2EuJifM = 0. 



2q 



N 



(105) 



Let us expand M and N in series of 



M = Mei + fME2 + ■■■ 

N = Nei + fNE2 + q^NEs + 



(106) 
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Substituting into ()104|) and equating successive powers of we obtain for the powers zero 
and two of g, respectively 

Mei = 0, (107) 



Me2 = [^ + ^i^J ^^^1- (108) 

Substituting ()106|) into ()105|1 we find that the equahty is identically satisfied at the zero 
order of q. For the orders and q^ we obtain, respectively: 

^^^^1 + -N'ei + Nei = 0, (109) 



Ne2 + In'e2 ~{} + ^) = - (^i' + I) Nei - (| + 2^^i) Me2. (HO) 

Let us note that the coefficient multiplying Ne3 is zero. 

The solution of ()109|1 bounded at infinity is the modified Bessel function of the first 
order Ki{s). Its expansion about s = is: 

A,,..c{i + |.,>|-£(l-2o04.„|-|;(5-2o)+...} (Ill) 

where C is an arbitrary parameter and 7 the Euler contant. 

From ()l()8j) we have Me2- Now the r.h.s. of Eq. ()11()|1 is known, since (jlllll is the so- 
lution of the homogeneous part of the equation (jllOj) . by the variation-of-constant formula 
we obtain the solution of (jllOj) : 

1 ln^(s/2) Ins /I ^- 
1 /9 7 ln2 ln2(2)\ 1 , , 



s V8 4 4 2 ^ 

From (fT06|l . (fToTjl . (fT08|l . (fTTT^ . and (fTT2|l . we have: 

^.c|l.£,£-i(l-2.)4.ni-£;(5-2.). 



2 



if!-2 + !!i!-!!!M|+...^. (114) 



s V8 4 4 
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6.3 Matching 

The interior solution for y ^ oo must match the exterior one for s ^ in an overlap 
domain. For the comparison we express the inner solution in term of s: 

^ f 1 1 , s 1 /ln2 1 , . ^ M , ^ 

M = ,3 - + -ln- + - ---ln, + A, + i?., +..., (115) 



iV = g<^- + -ln- + s Ing + Ai 

2 2 V 2 2 ^ 

ofl 1,9s Ins/, , , - 3 

+f |-- + - In^ - + - ^ln2 - Ing + 2A, + 2E., + - 

1 / ln2 2 In^g , / - 3\ \ 1 , , 

+- i — + ^~\nql2A, + 2Euj, + -\+A2j+---\. (116) 

A comparison of ()113p . ()114|) with pi5|) . ()116|) for small s shows that 

C = q , 2Ai + l = 7-ln-. (117) 

It is surprising that these two equalities are identical to the corresponding relations of 
Grant In fact, ()117p is satisfied at the order q and for the three first order of s in 
M and N. 

Then using (fTU^ . (fTUHIl and (fTTTj) we obtain 

Eo; = ^ fin- -7-0.115y (118) 



2 V q 

If we return to the physical variables, we have the dispersion relation: 

uj = (\n— --f -0.115] , (119) 



2E V qro 

that is, with 



g2 / 2 



where ro is the core radius ()59|1 and fo the velocity at r = tq of the relativistic classical 
fluid. 



7 Conclusion 

In the begining of this paper we consider the relations between the quantum nature of 
the condensate of an almost ideal relativistic Bose gas and its phenomenological properties 
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as a relativistic superfluid. The quantum origin of the description imposes inevitably the 
equation of state of the fluid (jlBj) which gives an acceptable velocity of sound. Incidently, 
we note that if we choose a potential with spontaneously broken symmetry (by changing 

into —m? in the velocity of sound would always be greater than 

The study of a cylindrical vortex has shown that this equation of state makes sense 
outside the core of the vortex whose the radius is approximately of the order of the healing 
length which indicates the onset of the quantum effects. 

Concerning the Kelvin waves, it is surprising that the mathematical calculations are 
so close the non relativistic ones. Yet the field equation differs from the Gross-Pitaevskii 
equation essentially in the second time derivative. Consequently the main equations (jHHI), 
fl89|) are different from the corresponding ones of The new dispersion formula ()120|) 
differs from the nonrelativistic one by the relativistic factor a/1 — where vq represents 
the velocity of the fluid in r = ro the radius of the core of the vortex. 

It is interesting to compare with the spinning global strings which behave like vortices 
in a superfluid medium. Ben-Ya'acov has given a covariant dispersion relation of the 
Kelvin waves. In the rest frame of the materiel medium, for an unit winding number and 
to < q, using (jH), it can be express as 



an helix at rest. E is not any more given by (j6(J|) but tends towards the particle mass in 
the non relativistic limit. 

The formula ()120|) and 1)1211) are different but, both are compatible at weak energy 
with the Grant formula (j21). However it is worth stressing that ()120|) is certainly the better 
dispersion relation for a weakly relativistic almost ideal Bose gas. The superfluid medium 
concerned by the dispersion relation ()121|1 is unusual as the velocity of sound is greater 



To conclude we note that relativity does not seem very useful in the laboratory exper- 
iments on the condensates but could be of interest in astrophysics. 
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(121) 



where Vi = - 



is the velocity along z of the Lorentz frame in which the pertubed vortex is 



than l/VS 
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